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Abstract 

The decays of the X(3872) to P-wave quarkonia are calculated under the assumption that it is a 
shallow bound state of neutral charmed mesons. The X(3872) is described using an effective theory 
of nonrelativistic D mesons and pions (X-EFT). We calculate X(3872) decays by first matching 
heavy hadron chiral perturbation theory (HHxPT) amplitudes for D*^ — > Xcj (tt", vttt) onto 
local operators in X-EFT, and then using these operators to calculate the X(3872) decays. This 
procedure reproduces the factorization theorems for X(3872) decays to conventional quarkonia 
previously derived using the operator product expansion. For single pion decays, we find nontrivial 
dependence on the pion energy from HHxPT diagrams with virtual D mesons. This nontrivial 
energy dependence can potentially modify heavy quark symmetry predictions for the relative sizes 
of decay rates. At leading order, decays to final states with two pions are dominated by the final 
state Xcitt'^tt'^, with a branching fraction just below that for the decay to Xci^r^. Decays to all 
other final states with two pions are highly suppressed. 
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The X(3872) is a novel charmonium state that hes very close to the D^D*^ 

threshold. The closeness of the X(3872) to this threshold (—0.6 ± 0.6 Me V) has prompted 



numerous authors to sue 



for reviews see Refs. 



?est that it is a C = 1 molecular bound state of D^D*'^ + D^D*^ 



5[. The strongest evidence for the molecular hypothesis is the large 
ratio of the branching fractions for the X(3872) decays to J/ip plus two or three pions 6| 

B.|X(3872) J/^.W] 

Br[X(3872) J/^tt+tt-] ^ ^ 

The final states have opposite G-parity which implies that the X(3872) does not have definite 
isospin. Another possibility, suggested in Ref. is that there are actually two nearly 
degenerate states with opposite G-parity. Measurements of the invariant mass distributions 
of the pions indicate that these decays proceed through J/ip P for the J/ipTc^ tt~ final state, 
and through J/ipu for the J/ipn^ 7r° final state, suggesting that the X(3872) couples to 
/ = and 1=1 channels with roughly equal strength. This is not possible for a conventional 
cc state. The observed branching ratio can be understood if the X(3872) is comprised 
primarily of neutral D mesons, which is expected since the charged D-D* threshold is far 
(8.5 MeV) above the X(3872) mass. 

If the X(3872) is a molecular state then the large distance behavior of the DD* wave- 
function is determined entirely by the binding energy and effective range theory (ERT) can 
be used to calculate some properties of the X(3872). ERT exploits the fact that at long 
distances the two-body wavefunction of the X(3872) must be 



e 



-■yr 



^x(r)(x , (2) 



where r is the separation of the meson constituents and 7 = ^/2BJlod*, where B is the 
binding energy of the X(3872), and udd* is D^-D*^ reduced mass. Parameter-free predic- 
tions for the decays X(3872) D^D^tt^ ^\ and X(3872) D^D^'y 9] have been obtained 
using ERT. Experimental measurement of these decays consistent with these calculations 
would confirm the molecular hypothesis. Unfortunately, data which can test these calcu- 
lations is not currently available. For example, the photon energy spectrum in the decay 
X(3872) D^D^'y near the kinematic endpoint can be computed in terms of the binding 
energy, but this spectrum has not been measured yet. The ratio of branching ratios flo| 



Br[X(3872) ^ L)°L)°7rO] 
Br [X (3872) ^ J/ 1^71+71-] ~ "'"-^'^ ' 



2+3.1 



(3) 



has been measured but cannot be calculated within ERT since the decay X(3872) 
J/ipTT^TT^ is sensitive to the short- distance structure of the X(3872). 

Many of the observed decays of the X(3872) are to final states that include conventional 
quarkonium such as X(3872) J/ip tt+tt-, X(3872) J/^ 7r+ tt" 7r°, and X(3872) ^ 
J/ip'y. These decays necessarily involve shorter distance scales where the description of 
the X(3872) as a loosely bound state of D and D* mesons is no longer valid. To make 
quantitative calculations for these decays one can use factorization theorems that separate 
long-distance physics {p ^ 7) from physics at shorter distance scales {p ^ 7) [ll|, Here 
7 can be determined from the binding energy B using the known masses of the D^, D*^, 
and X(3872). For m^o and m^*o, we use 1864.84 ± 0.18 MeV and 2006.97 ± 0.19MeV, 
respectively [ij]. For the X(3872) we use the mass obtained in measurements of decays to 
J/ipirir final states only, which yields mx = 3871.2 ± 0.5 MeV [l^. ^ Using these numbers 
we find B = 0.6 ± 0.6 MeV, so 7 ^ 48 MeV. The factorization theorems express the rate 
for X(3872) / (where / is the final state) in terms of the the rate for D^D*^ —>■ f times 
universal factors that depend on the binding energy of the X(3872). These factorization 



theorems have been used to analyze X(3872) decays into J/ip plus light hadrons in Ref. |15 |. 

Recently, Ref. [l^ has studied the decays X(3872) — > Xcj(7r°, tttt) and argued that mea- 
surement of the partial rates for these decays can discriminate between conventional char- 
monium and molecular interpretations of the X(3872). Decays to Xcj are particularly in- 
teresting since the Xcj are in a heavy quark multiplet and therefore heavy quark symmetry 
can be used to make predictions for the relative sizes of the partial rates. 

In this paper, we study the strong decays X(3872) — > y^jfTT^, tttt) using a combination 
of heavy hadron chiral perturbation theory (HHxPT) |l7, 18, 19| and X-EFT, an effective 



field theory for the X(3872) developed in Ref. 2^. HHxPT is an effective theory for 
low energy hadronic interactions of mesons containing heavy quarks that incorporates the 
constraints of the heavy quark and chiral symmetries of QCD. We use HHxPT to calculate 
the D^D*° — > Xcj(7r°, tttt)^ transition amplitudes. X-EF' 
developed for shallow nuclear bound states in Ref. [21 . 



20 [generalizes the effective theory 



HQ 



231]. In this theory, modes have 



^ Interpretation of mass measurements in decays to tt" final states is complicated by the possibility of 

two distinct states as well as threshold enhancements due to the large scattering length in D'^D*^ 

channel. We thank E. Braaten for discussion on this last point. 

^ Here and below decays related by charge conjugation are implied. 
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momenta of order 7 so the degrees of freedom are non-relativistic charmed mesons and non- 
relativistic pions. The decays X(3872) — > Xcj(7r°, tttt) proceed through local operators that 
couple the and D*^ to the Xcj and one or two pions. We use the HH^PT calculation 
of the D^D*^ Xcj(^°;7r7r) amplitudes to fix the coefficients of these operators in the X- 
EFT Lagrangian, then calculate the decays of the X(3872). The resulting decay rates are 
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121. 



consistent with the factorization theorems for X(3872) decays proven in Refs. 

In the leading order (LO) diagrams for D^D*^ Xcjtt^, virtual charmed meson prop- 
agators introduce dependence on the pion momentum which modifies the usual de- 
pendence expected for a P-wave single pion decay. This can give significant correc- 



tions to the relative rates predicted in Ref. 16|. We also calculate the differential decay 



rates for X(3872) Xcjmr. At LO in HHxPT, amplitudes for ^ Xcjtt" and 

D^D*^ —>■ Xcj TT TT are given in terms of a single unknown coupling in the HHxPT lagrangian. 
Therefore, at LO, the normalization of the X(3872) XcJ^t^^ decay rates can be given in 
terms of the rates for X(3872) — > Xcj tt". We find that X(3872) — xa tt" 7r° is the dominant 
decay mode with two pions in the final state, with a branching fraction that is just 0.6 times 
the branching fraction for X(3872) Xa 7r°. All other decays to final states with two pions 
are considerably smaller. In some of the tree level diagrams for D^D*^ Xcitt'^tt^, a virtual 
meson can go nearly on-shell, with the potential divergence cutoff by the width of the 
£)*o rjij-^g resulting contribution for the decay rate is enhanced by ~ 27r£^^/r£i.o ~ lO'*, 



state (about 200 MeV) and r^*o, 
24l |. This enhancement only occurs 



where E^^ is the typical energy of the pion in the fina 
the width of the Z)*°, is estimated to be about 68 keV 
for the Xci7r°7r° final state, total partial decay rates to other final states with two pions are 
smaller by orders of magnitude. 

The LO HHxPT lagrangian for mesons containing heavy quarks or antiquarks at rest is 

-f i:i[Hl{tDo)abH^] + 9 "^[Hl a ■ Aai^H,] + ^Ti[HI a' a'] . (4) 



We use the two component notation of Ref. 



2M[. The field Ha is given by 



Ha = Va-a + Pa, (5) 

where Va annihilates D* mesons and Pa annihilates Da mesons. The subscript a is an SU{3) 
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index, and a = 1 for neutral D mesons. The field for antimesons is 



Ha = -Va-a + Pa. (6) 

The field Aab is the axial current of chiral perturbation theory, Aab = —'Vnab/fn + where 
Jtt is the pion decay constant and Hah are the Goldstone boson fields. The axial coupling, g, 
is taken to be 0.6 and Ah is the hyperfine splitting. 

To calculate the amplitude for D° D*^ — Xcj ti"" we need to include the Xcj in the HH^PT 
lagrangian. The Xcj are degenerate in the heavy quark limit and form a heavy quark 
multiplet. The fields for Xcj have been constructed in a covariant formalism in Ref. [25|. In 
our two component notation, the Xcj fields are represented by 

X' = CT^ X'^ 

= (x^ + 4 ^''"xl + ^Xo) . (7) 



V2 v^' 

The transformation properties of the heavy quark and quarkonium fields under the various 
symmetries are: 

rotations: Ha^UHaW Ha^UHaW x" ^ R'^U^X^U 

heavy quark spin : Ha SHa Ha HaS"^ X* ~^ Sx'^S'^ 

parity : Ha -Ha Ha ^ -Ha x' X' 

charge conjugation : if„ ^ asi^Jas Ha^ a2Hja2 x' ~^ -(^2{x'Y(^2 = X' 

SUl{3) X SUr{3) : Ha^HbVl Ha ^ VatHb x'^X' 

Note that under charge conjugation, the components of Ha transform as Va V a and 
Pa Pa- Therefore, the C even combination of D mesons that couples to the X(3872) 
is ViPi + ViP. In the first line, are 2 x 2 (3 x 3) rotation matrices related by 

U'^a'^U = R'^aK In the second line, 5* is a rotation matrix acting only on the heavy-quark 
spin, 5 is a rotation matrix acting only on the heavy-antiquark spin and in the last line Vab 
is an SU (3) matrix. 

The lagrangian coupling the Xcj to heavy mesons is 

= t^Trix^'Haa' Ha] + jTiix^'Haa^ Hb^j^A'ab + h.c. . (8) 

The leading contribution to the amplitude for D*^ Xcj comes from the tree level 
diagrams in Figs. [1^), b), and c). In these diagrams, the coupling of XcJ to D mesons comes 




FIG. 1: Feynman diagrams contributing to the D D XcJT^ amplitude. The sohd Une is a D 
meson, the double line is a D*^ meson, the dotted line is a vr" meson, and the thick solid line is 
the Xc meson. 

from the first operator in Eq. ([8]). Fig. [1^) contributes to Xco decays, Fig. [T)d) contributes to 
both Xco and Xc2 decays, and Fig. [T]^) contributes to Xci decays only. In the HHxPT power 
counting, one takes Ej, ~ ~ A// ~ Q, and expands amplitudes in powers of Q. The 
graphs in Figs. [T^)-c) are 0{Q^) since the coupling of the D mesons to Xcj has no derivatives 
and is therefore 0{Q^)^ the axial coupling of the D mesons is 0(Q), and the heavy meson 
propagator is 0{Q^^). Fig. [TJi) shows the contribution to the amplitude from the second 
operator in Eq. ([8]). This contribution is 0{Q) because this operator contains one derivative. 
Chiral loop corrections to the LO diagrams in Eq. ([1]) are suppressed by and therefore 
subleading to the diagram in Fig. [T]i). It is also clear that higher dimension operators that 
couple the D mesons to Xcj either contain two spatial derivatives or an insertion of the 
light quark mass matrix and are suppressed relative to the leading operator in Eq. ([H]) by 
at least . Therefore the diagrams in Fig. [1] give the complete HHxPT calculation to 
next-to-leading (NLO) order. 
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We obtain 



2 CD ■ Ptt 



V2 



9 9i 



+ 



+ 



V2 



9 91 



1 



V2. 



(9) 



The results for the charge conjugate initial states are identical. These amplitudes must be 
multiplied by a factor of ^rno tud* f^^x^j to account for the non-relativistic normalization of 
the HHxPT fields. 

The virtual charmed mesons in Figs. [T^)-c) are off-shell by ~ -E^, where E^^ varies from 
432 MeV to 305 MeV in the three decays. X-EFT contains non-relativistic , D*^, and 
7r°, and total energies of these particles are assumed to be within a few MeV of the D^D*^ 
threshold. Thus we must match the tree level diagrams in Fig. [1] onto local operators in 
X-EFT. We will illustrate this for the decay to Xco? the generalization to other decays is 
straightforward. The operator in X-EFT which reproduces the amplitudes for Z)*° 
Xco tt" and D° 



Xco 7r° is 



.C*x,o(-E'7r) 2 



AcO : 



(10) 



where we have dropped the SU (3) labels on the D meson fields that are implicitly neutral. 
Note, when evaluating the amplitude we must multiply the naive Feynman rule obtained 



from Eq. (ITO!) by 1^16 ^xco "^d f^D* to account for the normalizations of the fields in the 
effective theory, which are all nonrelativistic, including the pion. The factor of 1/ (4y^m7) 
is included to compensate for the differing normalizations in the effective field theory and 
HH^PT. Comparing with Eq. (times y/iriD mo* ^xco); C'x,o(-^7r)) 



C^,o{E^) = V2gg, 



+ 



+ 



Cl 



(11) 



AE^-Ah ' AE^ + AhJ ' V2 
Now we consider the decay X(3872) — > Xcjt^^ mediated by this operator. The one- 
loop diagram contributing to this decay is depicted in Fig. O The decay matrix element 
is obtained from the diagram in Fig. [2] after dividing by the wavefunction renormalization 
factor jiDD* I V27i"7 [20]. After summing over both channels, the result is 

1 



A<(X(3872)^Xcovr' 



— r V"^xO C^fiiE^) ^Pn-ex. (12) 




FIG. 2: Feynman diagram contributing to the X(3872) XcJT^^ decay amplitude. The circle with 
cross represents the interpolating field for the X(3872). 



In evaluating the diagram in Fig. [2] we have used 

d^l 1 



2/z 



DD* 



A^DD* 



(A 



PDS 



7) 



It is straightforward to evaluate the rates: 

nn _ 27r7 



r[X(3872) ^ Xcjtt' 



m 



XcJ 



PDS) 



/i2,^.C2(ApDs) mx 727r/2 
where ao = 4, ai = 3, 02 = 5, and the functions C^^j{Et,) are given by 

Ci 



(13) 



(14) 



^1 



Using these results we can calculate the partial rates for r[X(3872) 
(ci = 0) we find (denoting r[X(3872) ^ Xcjtt^] = V j) 



(15) 

XcJ A- At LO 
(16) 



Fo : Fi : F2 :: 4.76 : 1.57 : 1. 

The energy dependence from the virtual pion propagators makes a significant change in the 
predictions for the relative sizes of the partial rates. The NLO predictions depend on the 
ratio Ci/gi which is undetermined. The ratio ci/gi has dimensions of inverse mass, when 
ci/gi is varied from (100 MeV)-^ to (1000MeV)~\ we find: 

Fo : Fi : F2 :: 3.01 : 1.06 : 1 Ci/gi = (100 MeV)"^ 

Fo : Fi : F2 :: 3.49 : 1.20 : 1 Ci/gi = (300 MeV)"^ 

Fo : Fi : F2 :: 3.76 : 1.28 : 1 Ci/gi = (500 MeV)"^ 

Fo : Fi : F2 :: 4.11 : 1.38 : 1 Ci/gi = (1000 MeV)"^ . (17) 
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a) 



b) 



0,2 



c) 



1 d) 



0,2 



e) 



f) 







g) 



FIG. 3: Feynman diagrams contributing to the D^D*^ XcJ^r'^vr'^ amplitude. Diagrams with the 
final state pions crossed are not shown. The numbers denote the angular momentum quantum 
number of the Xcj produced in the decay. 

For the largest values of ci/gi, we obtain predictions for the relative sizes of the partial 
rates similar to those found in Ref. [igI. As Ci/gi decreases, the predictions tend to those 
appearing in Eq. ( flGll . Experimental measurement of the relative sizes of the partial rates 
for decays to Yc,j7r° can be used to determine ci/gi. 



Refs. 
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12| derived factorization theorems for X(3872) decaying into quarkonium states 



and light hadrons. Note that the first term in the decay rates in Eq. (fT4l) can be written as 

I (0| ;7f ^«(^) P + V'P) |X(3872, X))\\ 



/^i)_D*^o(^PDs) 



(18) 



where e{X) is a polarization vector. Thus we see that the expressions for the partial widths 
factor into the product of long-distance matrix elements, given in Eq. (ITSj) . times short- 
distance coefficients, which are proportional to the cross sections for D*^ Xcjt^^- This 



12l |. A main point of this paper is 



is the content of the factorization theorems of Refs. [11 
that HHxPT can be used to calculate the short- distance coefficients in these factorization 
theorems. 

Next we turn to the calculation of the three-body decays, X(3872) Xcj ^r*^ vr'^. In Fig. [3] 
we show the diagrams contributing to the transition amplitudes for Z)*° Xcj 7r°. 
These amplitudes are matched onto local operators in X-EFT, then these local operators 
are used to calculate the decays X(3872) Xcj t^^, in a manner identical to the procedure 
used to calculate X(3872) — > Xcj 7r°. Here we give our results for the differential partial decay 
rates. 

The differential decay rate for X(3872) Xco'^^'^^ is given by 
dEidE2 /il)z)Xo(ApDs) 18 

where 

1 



F{Ei, E2) = {E2 — E 



[(El -t- ^h){E2 + Ah) (El + E2 + Ah) 

3 , Ah 



(20) 



El E2 (El + E2- Ah) El E2 (El + Ah) (E2 + Ah) 

Here, Ei{pi) refers to the energy (three- momentum) of one of the 7r° and Pi = \pi\. The 
partial decay rates are symmetric under 1^2. In these decays Ah is equal to the neutral 
hyperfine splitting. Ah = Ah = Tnij*o — m/jo = 142 MeV. 

The rates for X(3872) Xcj'^~^'^~ can also be calculated. The differential rates for 
X(3872) — i> XcjTT^TT^ are bigger by a factor of two than the corresponding rates for 
X(3872) —>■ XcjTT^Tr" due to isospin factors in the diagram. Furthermore, the factor of 
Ah can be A^, A^+ = mD*o - mD+ = 137 MeV, or Ajf = mc,*+ - ttido = 145 MeV, de- 
pending on the diagram, so the expressions above (times 2) can only be used if one neglects 
isospin breaking in the heavy meson hyperfine splittings. 

The differential decay rate for X(3872) — > Xcitt^tt^ obtained from diagrams in 
Figs. [3)3), c),e), and f) is 

d^T[X -^Xcin^Tv^] _ 27r7 vr (7 ^X.i 



dEidE2 filo'CoiApDsy 3 

X pI pI {Bf^ + Bl^) + (pi • p2f (3 ^2 + 2 A (E12 + E21) + 2 E12 E21) 



(21) 
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where 

Bi2 = B{Ei,E2 



El E2 {El + E2){Ei + A^^) {E2 + A^^) ' 
2 2Ei + Ah 



{El + Ah) {E2 - Ah) Ei {Ei + E2) (^i + Ah) ' 
B21 = B{E2,Ei). (22) 

In Figs. [3b) and[3]F), the virtual meson can go on-shell leading to the poles of Ei — Ah 
in the expression for the decay rate. In the case of vr^ tt" final states, the virtual charged 
D meson cannot go on-shell but E^ — Ah can be as small as 2 MeV. The pole in the 
decay amplitude is regulated by putting the external D*^ on its complex energy mass shell, 
E£,*o = AH — iT£,*o/2. ^ The width of the D*^ has not been measured. If one uses HHxPT to 
relate the width of the D*^ to the measured width of the D*~^ and the measured strong decay 

n 

branching fractions for both D*^ and D*^, one finds To*o = 0.068 MeV [241]. Since this scale 
is orders of magnitude smaller than all other scales in the problem, we will keep only terms 
which do not vanish in the limit F^i.o 0. For terms with single poles in E^ — Ah, this 
amounts to using a principal value prescription. When integrating over the double poles, we 
replace {Ei — Ah)"^ — * {Ei — Ah)"^ + F^,o/4 in the denominators of Eq. (!2T1) . then perform 
the energy integrals using the formula 

where the elipsis denotes terms that can be neglected when F <^ xo,xi. Keeping only the 
first term in Eq. (!23l) is called the narrow width approximation which turns out to be an 
excellent approximation for X(3872) — > Xci^t^tt^. 

One might also worry about contributions from interference terms that contain single 
poles in Ei and E2. These are of the form 

Re (1^ W^—T^l^ W^—^) (24) 



El- Ah- tT/2 E2-AH + 

{Ei-Ah){E2-Ah) , FV4 



((Ei-A^)2 + FV4)((E2-Ah) + FV4) {{Ei-Ah)^ + T^IA){{E2-Ah) + T^IA) 



We thank Eric Braatcn for suggesting this. 
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The first term on the righthand side clearly yields two single poles which should be dealt 
with via the principal value prescription in the limit Two — > 0. The second term only has 
a nonvanishing contribution (in the limit F^.o — > 0) with support at E1—E2 — A/^, which 
lies outside available phase space. 

Finally, the differential decay rate for X(3872) — >• Xc27i"°7i"° is 



dEidE2 /x|,^.Co(ApDs)' 3 

X pIpI (C2 + C(L>i2-L'2l)) + ^PL + ^2l) + ^A2i?2l) , 

• P2f (-^ + C {D^2 - D2i)) - ^{Dl, + Dl,) + ^ /^2i) 



(25) 



where 



C = C{Ei,E2) 



(^1 + Ah) (^2 + Ah)(^i + E2 + Ah) 

^-"^^^-'^-^^ (^./a.)^. ' 

D21 = D(E2,E^). (26) 

The normalization of the three-body decay rates is unknown but can be given in terms 
of the LO chiral perturbation theory expressions for the two-body decay rates computed 
earlier. Integrating over three-body phase space we find 

•Br[X(3872)^Xco7r0 7rO]\ ^g-^^^-,^ 

/ LO 

/Br[X(3872) ^ xa vr° vr"] \ ^ ^ ^ . 



01 / =6.110- 



Br[X(3872)^Xco7rO 
lr[X(3872) - 
V Br[X(3872) 
/ Br[X(3872)^Xc2vr"7r°] \ _ , 

[ Br[X(3872)^Xc2vrO] J.^" ' ^ ^ 

Here we have used mx = 3871.4 MeV and = ^Trfj^ = 1659 MeV. The subscript LO 
emphasizes that this result only holds at LO in xPT. At NLO the prediction depends on 
the unknown parameter, ci. 

We have not numerically computed the corresponding branching fraction ratios for 
X(3872) XcJ'^'^'^~- For final states with J = or 2, isospin violation is a small ef- 
fect so we expect 

/ Br[X(3872)-^Xco7r+7r-] \ ^ ^ /^ Br[X(3872) ^ Xeo tt" vr"] \ ^ ^ , 

V Br[X(3872)^Xco7rO] J \ Br[X(3872) ^ Xco vr"] ' 

/ Br[X(3872)^Xc2 7r+7r-] \ ^ / Br[X(3872) ^ Xc2 vr" 7r°] \ ^ 

1, Br[X(3872)^Xc2vrO] ^ % Br[X(3872) ^ Xc2 vrO] ~ ^"^ " ^^'^ 
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For final states with Xci tt"*" vr~, isospin breaking effects are crucial since the virtual D mesons 
cannot go on-shell and therefore the first term in Eq. fl23|) which is responsible for the 
enhancement of X(3872) Xci tt*^ is absent. However, the decay Xci tt^ is enhanced 
relative to Xco,2 tt^ tt^ decays by factors of Et,/ {E^, — A^) ~ 10^, so we expect 

/ Br[X(3872)^Xcovr+7r-] \ ^ 3 
^ Br[X(3872)^Xco^°] Ao ^ 

(This is the same order of magnitude as the subleading terms in Eq. fl25]) for the decay 
X(3872) — > Xci 7r°7r°.) In the absence of available data on X(3872) — > XcJ^^t^, we will not 
attempt a more precise estimate in this paper. 

Except for the decay, X(3872) —>■ Xcitt^tt", the decays to final states with Xcj plus two 
pions are quite small compared to the two-body decays. These decays are suppressed by 
factors of g^p^ / as well as limited phase space. In the case of X(3872) Xcitt^tt", 
there is an enhancement by a factor of 2TiET^/T £,,0 ~ 10^, which significantly enhances the 
decay rate, making it comparable to r[X(3872) — > Xci tt^]. This enhancement is due to the 
virtual D meson propagators in Figs. [3] c) and f), which can go on-shell with -k^ -k^ final 
states. The potential divergence in this graph is regulated by the tiny width of the D*^. 
Ref. [16|] observed that the three-body decays to xa tt vt should dominate over decays to final 
states with Xco or Xc2 when the X(3872) is modeled as a ^Pi charmonium, however, in this 
scenario the predictions for the ratios in Eq. fll6m7p are quite different. For the hadronic 
molecule scenario, Ref. [16] applies the QCD multipole expansion and finds that the rate for 
X(3872) — i> XcoTTTT vanishes, but does not make quantitative predictions for the three-body 
decays involving Xci and Xc2- 

In summary, we have calculated the decays of the X(3872) to P-wave charmonium us- 
ing a combination of X-EFT, which treats the X(3872) as a shallow bound state of the 
nonrelativistic D mesons, and HHxPT. Our methods reproduce established factorization 
theorems for decays of the X(3872). HHxPT can be used to calculate the short-distance 
coefficients in these factorization theorems. For two-body decays to Xcj tt" we noted that 
the p\ dependence that is usually assumed for P-wave single pion decays can be modified 
by the virtual D meson propagators in the diagrams shown in Fig. [2l This can affect heavy 
quark symmetry predictions for the relative sizes of the partial decay rates. Finally, we find 
that at LO the partial rates for the X(3872) Xci^t^tt^ is only slightly smaller than the 
rate for X(3872) Xci^r^, while other decays to final states with XcJ and two pions are 
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orders of magnitude smaller. 
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